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Hydrodynamics of Cloud Collisions in 2D: 

The Fate of Clouds in a Multi-phase Medium^° 

Francesco Miniati T. W. Jones^’^, Andrea Ferrara and Dongsu Ryu 

ABSTRACT 

We have studied head-on collisions between equal-mass, mildly supersonic HI clouds 
(Mach number 1.5 with respect to the background medium) through high resolution 
numerical simulations in two dimensions. We explore the role of various factors, 
including the radiative cooling parameter, r] = Trad/ TcoII ( TcoIi=Rc/vc), evolutionary 
modifications on the cloud structure and the symmetry of the problem. Self-gravity is 
not included. Radiative losses are taken into account explicitly and not approximated 
with an isothermal adiabatic index 7 ~ 1 , which, in fact, leads to very different results. 
We assume a standard two-phase interstellar medium (ISM) model where clouds are 
characterized by a temperature Tc = 74iF and number density ric = 22cm~^ and are 
in pressure equilibrium with the surrounding warm intercloud medium (WIM), with a 
density contrast y = pdPi = 100. In particular, we study collisions for the adiabatic 
(7 ^ 1) and radiative (7 = 0.38) cases which may correspond to small {Rc < 0.4pc for 
assumed WIM) or large {Rc ~ l.bpc) clouds, respectively. In addition to a standard 
case of identical “non-evolved” clouds, we also consider the collision of identical clouds, 
“evolved” through independent motion within the intercloud gas, over one crushing 
time before collision. This turns out to be about the mean collision time for such 
clouds in the ISM. The presence of bow shocks and ram pressure from material in the 
cloud wake, significantly alters these interactions with respect to the standard case. In 
some cases, we removed the mirror symmetry from the problem by colliding initially 
identical clouds “evolved” to different ages before impact. In those cases the colliding 

^Department of Astronomy, University of Minnesota, Minneapolis, MN 55455 
^Dipartimento di Astronomia, Universita of Firenze, 50125 Firenze, Italy 
^Osservatorio Astrofisico di Arcetri, 50125 Firenze, Italy 

^Department of Astronomy & Space Science, Chungnam National University, Daejeon 305-764, Korea 

^Department of Astronomy, University of Washington, Box 351580, Seattle, WA 98195-1580 

®e-mail: min@msi.umn.edu 

^e-mail: twj@astro.spa.umn.edu 

®e-mail: ferrara@arcetri.astro.it 

®e-mail: ryu@sirius.chungnam.ac.kr 

^°Published in The Astrophysical Journal, December 10, 1997 



- 2 - 


clouds have different density and velocity structures, so that they provide a first insight 
on the behavior of more complex interactions. 

In our adiabatic collisions the clouds are generally disrupted and convert their gas 
into the warm phase of the ISM. Although the details depend on the initial conditions, 
the two colliding clonds are converted into a few low density contrast (y ~ 5) clumps 
at the end of the simulations. 

By contrast, for symmetric radiative cases we find that the two clouds coalesce, 
and there are good chances for a new massive cloud to be formed. Almost all the 
initial kinetic energy of the two clouds is radiated away during such collisions. On 
the other hand, for both adiabatic and radiative collisions, symmetry breaking leads 
to major differences. Most importantly, asymmetric collisions have a much greater 
tendency to disrupt the two clouds. Portions of individual clouds may be sheared 
away, and instabilities along the interfaces between the clonds and with the interclond 
medium are enhanced. In addition, radiative cooling is less efficient in our asymmetric 
interactions, so that those parts of the clouds that initially seem to merge are more 
likely to re-expand and fade into the warm intercloud medium. Since the majority of 
real cloud collisions should be asymmetric for one reason or another, we conclude that 
most gasdynamical diffuse cloud collisions will be disruptive, at least in the absence of 
significant self-gravity or of a significant magnetic field. 


Subject headings: ISM: kinematics and dynamics - hydrodynamics 



- 3 - 


1. Introduction 

Interstellar cloud collisions (CCs hereafter) are important for the dynamical evolution of 
galaxies. Indeed, CCs turn out to be relevant in a large variety of important processes occurring 
in the interstellar medium (ISM), such as star formation, dissipation of kinetic energy and gas 
phase transitions. In addition, they affect the gaseous structure and the energy budget of galaxies, 
along with the mass spectrum of diffuse ISM clouds and their evolution. 

CCs were first studied by Stone (1970a, 1970b), who was mainly concerned with the evolution 
of interstellar clouds and especially of their structure. In his pioneering papers he found that, 
despite the very high compression undergone by clouds during CCs, star formation is not 
enhanced. Instead, he found that colliding clouds lose mass and that large scale perturbations 
and internal motions (approximately radial) form within the clouds, which do not decay before 
the next collision. Thus, cloud hydrostatic equilibrium is severely compromised, suggesting the 
idea of clouds with smooth density distribution profiles instead of (clouds as) sharply bounded 
objects (Stone 1970b). Smith (1980) performed ID numerical simulations and concluded that, 
in a dusty ISM, low Mach number (< 5) collisions are the most likely to trigger star formation. 
Further investigations have been carried out by Gilden (1984) and Lattanzio et al. (1985, LMPS 
hereafter) using 2D and 3D models, respectively. These authors studied the role of various 
parameters involved in the collisions, such as the relative size of the clouds, their masses and 
impact parameter. They generally concluded that CCs more often lead to disruption than to 
coalescence or gravitational collapse of the clouds. This again has been interpreted as evidence 
that the ISM should be described as a continuous distribution of gas streams and turbulent eddies 
(LMPS, Hunter et al. , 1985). That contrasts with most dynamical models for the ISM that depict 
it more simply as pressure-confined neutral cold clouds (CNM) moving through a warm diffuse 
inter-cloud medium (WIM) (see e.g., Norman and Ferrara 1996, Vazquez-Semadeni et al. 1995 
and references therein). This point is so important that further investigation of the final fate of 
the clouds is worthwhile. In fact, it turns out that CCs are very frequent, with a rate of about 
once every 10® — lO^yr for individual clouds, leading to a total of ~ 1 cloud collision every 100 yr 
in the Galaxy (Stone 1970, Spitzer 1978, Hausman 1981). Their outcome could profoundly affect 
the mass interchange of the various phases that are believed to exist in the ISM. 

CCs may also be responsible for the buildup of the observed mass spectrum of diffuse 
clouds (Dickey & Garwood 1989, Solomon & Rivolo 1989, Oort 1954, Field & Saslaw 1965, Field 
& Hutchings 1968, Pension et al. 1969, Gowie 1980, Hausman 1982, Humphrey & Scalo 1983, 
Struck-Marcell & Scalo 1984, Fleck 1996, Mousumi & Chanda 1996). Early models suggested 
a hierarchical scenario in which small clouds, formed in HII regions out of much larger clouds 
through star formation, are supposed to undergo repeated inelastic collisions and coalescence, 
thus engendering new clouds of larger mass (Oort 1954, Field &: Saslaw 1965). Yet, an accurate 
model should take into account the detailed outcome of CC processes. In their model Field &: 
Saslaw (1965) assume, for example, that each CC leads to coalescence, an uncertain assumption as 
discussed above and by several authors including Stone (1970a, 1970b), Hausman (1981), Gilden 
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(1984) and Klein, McKee & Wood (1994, KMW hereafter). 

In addition, CCs are relevant for the energy budget and gaseous structure of a galaxy. Clouds 
are accelerated in HII regions (Field & Saslaw 1965) and by supernovae. The turbulence pressure 
associated with their motion determines the vertical structure of the cold neutral phase in the 
ISM (McKee 1990, Ferrara 1993, Norman &: Ferrara 1996). The bulk motions, on the other hand, 
are dissipated through inelastic collisions (Spitzer 1978). It is also worthwhile to mention that 
the amount of energy dissipated during CCs depends on the dust content and on the metallicity 
of the ISM, as discussed by Ricotti, Ferrara & Miniati (1997, RFM hereafter), who determined 
the dependence of the elasticity of a collision on parameters such as relative cloud velocity, size, 
metallicity and dust-to-gas ratio. Their paper gives a useful overview of the characteristics of 
collisions and represents a complementary study to the present one. 

The present paper is devoted to a further investigation of the consequences of CCs. In 
particular, we hope to determine the fate (i.e. survival, dispersal in the WIM, coalescence or 
fragmentation) of clouds in a two-phase ISM with canonical properties. In addition, we investigate 
the evolution of the cloud kinetic energy. This could be retained by the clouds as thermal energy, 
radiated away during the collision, or transferred to the WIM. Very recently some attention has 
been devoted to the hydrodynamical details that characterize the flow structure in CCs (KMW, 
Kimura & Tosa 1996). KMW have pointed out important instabilities that develop in the collision 
process. In particular, they have shown that small amplitude surface perturbations may lead to 
the development of the “bending mode instability” (Vishniac 1994) in the cloud shocks, which, in 
turn, causes the colliding material to form filamentary structures. Our work extends the KMW 
study by adding individual cloud evolution through the intercloud medium, before the collision. 
This increases the realism of the simulations by taking into account the role of cloud bow shocks 
and wakes and of nonuniform cloud structures. By colliding two clouds that have undergone 
different amounts of prior evolution, we introduce a simple means to relax the symmetry imposed 
in most previous calculations. This turns out to alter dramatically the result of the interactions 
in ways that extend beyond the direct influence of the bending mode. We still limit for now our 
attention to head-on 2D gasdynamical collisions between clouds of (initially) equal mass and Mach 
number. 

The plan of this paper is as follows. In §2 we describe the physical problem under investigation 
and discuss the assumptions made. In §3 we describe the computational setup; §4 presents the 
main results of this study; discussion is given in §5 whereas summary and conclusions are in §6. 


2. Description of the problem 

It is known {e.g., , Stone 1970a, 1970b, KMW) that collisions between identical homogeneous 
clouds are generally characterized by three main evolutionary stages as shown in Fig. 1. These 
are (1) a compression phase in which a shock wave forms in each cloud and propagates through 
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it, converting cloud kinetic energy into thermal energy. If radiative cooling is efficient, a fair 
fraction of this thermal energy is radiated away. Otherwise, the pressure becomes high enough to 
drive shocks outward from the center of the collision. At the end of the compression phase the 
gas is highly compressed into a thin (compared to the initial cloud size) layer, whose thickness 
depends on the amount of radiative losses (Fig. Ha and top panel of Fig. |a). (2) A re-expansion 
phase begins when shocks generated during the collision emerge into the inter-cloud medium. 
During this phase a rarefaction wave propagates backwards into the clouds, forming a central low 
pressure and low density region (Fig. lb and 2). Also, a fast sheet of ejecta emerges from the 
contact discontinuity between the merged clouds. This structure, which appears like a “jet” in 2D 
collisions with mirror symmetry between the clouds, is subject to Kelvin-Helmholtz instabilities 
(KHIs) (Fig. Also, during the re-expansion phase, a reverse shock forms and interacts with 
the expanding cloud material. As will be shown in §^, the qualitative and quantitative details of 
this phase strongly depend on the adiabaticity and symmetry of the collision. (3) There follows a 
collapse phase when the expansion is halted by the external surrounding medium. The pressure 
inside the clouds is by now much lower than in the external ambient medium, so that the cloud 
material is driven back toward the merged cloud core. As pointed out in the literature, during 
this phase Rayleigh-Taylor instabilities (RTIs) become very important, especially in disrupting 
the cloud surface and ablating material away from it (Fig. ||c). The details of the evolution of 
these properties are obviously closely related to the previous, re-expansion phase and thus depend 
as well on the adiabaticity of the collision. Under some circumstances it appears that one should 
add (4) a dispersal phase, since the original clouds may be largely converted to WIM (Fig. |d). 

In the absence of self-gravity and for Mach number high enough that Mach scaling applies 
(KMW) the primary parameter controlling head-on collisions between identical clouds is 


^rad 
^c-^c 


( 2 - 1 ) 


where Nrad is the radiative cooling column density, while Uc and Rc are the cloud number density 
and radius, respectively. Note that this definition for r/ agrees with the one in eq. (5) of RFM; 
their adjustable parameter a is here taken equal to unity. Thus when comparing the two sets of 
results our r\ should be divided by a~^ ~ 3. If we express in terms of the radiative cooling 
time Trad, Uc and the individual cloud velocity Vc (A^rad = ncVcTrad)-, 


'^rad ^rad 

^ = ~5~T~ =-’ 

where, according to Spitzer (1978), 


3 kT 
2 nA(n, T, Z) 


( 2 - 2 ) 


(2-3) 


and A(n, T, Z) is the interstellar cooling function depending on number density n, temperature 
T, and metallicity Z. We have introduced the collision time TcoU = Rc/vc, which is a natural 
timescale for CCs and is about the time over which the compression phase occurs. The collision 
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is adiabatic if 77 ;^> 1, radiative if r/ ~ 1 and isothermal if ?? -C 1. Since according to eq. Trad 
only depends on the density and temperature of the cloud (and on the metallicity, Z), then, from 
eq. 2-2 rj oc VcjRc and we can infer that small, supersonic clouds undergo primarily adiabatic 
collisions, whereas interactions between large, slow clouds are mostly radiative or isothermal. 
Since in the adiabatic case TcoU<^ Trad, there is not enough time for the thermal energy to be 
radiated away during the collision. On the other hand, in the radiative regime Trad ~ TcoU, and 
during the collision a significant fraction of the energy associated with the clouds is converted into 
radiation. This fraction is even larger in the isothermal case. 


Occasionally in calculations such as the present ones, strong radiative cooling is taken into 
account approximately by setting the gas adiabatic index 7 ~ 1. That is, the flow is assumed 
to be inherently isothermal. This method allows much greater compression in the flow than the 
usual 7 = 5/3 case, an effect similar to that expected during a radiative compression. However, 
this approach is not appropriate for the problem at hand. As long as total energy is conserved, 
even a 7 ~ 1 gas will have a substantially higher pressure behind a strong shock than in front 
of it. That leads in the present situation to strong pressure gradients that drive gas flows out 
of the impact region during a collision. On the other hand, when radiative cooling is included 
properly, the thermal energy is removed from the interaction region (during the collision) before 
the re-expansion takes place, significantly reducing forces that drive gas away from that region. 
We carried out comparison simulations using 7 ~ 1 to approximate strong radiative cooling. They 
had properties very different from the properly radiatively cooled flows, actually resembling more 
nearly the adiabatic flows for the reasons already mentioned. 

Considering weakly radiative or fully adiabatic cases, on the other hand, it is worthwhile 
to notice that we performed several tests with 1.5 < p < 00 and did not find any substantial 
differences among them either in qualitative or quantitative terms. This means that even though 
our adiabatic simulations are characterized by 77 S> 1, they also represent reasonably well cases 
with 77 > 1.5. 

In this paper we restrict our study to head-on collision of neutral hydrogen (HI) clouds, 
and only consider supersonic clouds. In fact, according to Spitzer (1978), supersonic clouds 
should be the most common case for the ISM. We neglect for now the magnetic field. KMW and 
RFM included, in their calculations, a magnetic pressure term corresponding to an initially weak 
(B r\j l^G) magnetic field in order to limit the extraordinary compression otherwise occurring 
during radiative or isothermal collisions. However, a full MHD simulation was lacking. In complex 
flows, shear is typically more important to magnetic field behavior than compression (e.g., Frank 
et al. 1996, Jones et al. 1996), so full MHD may be expected to behave differently. We ignore any 
thermal conduction effects, as well as self-gravity in the present simulations. 

We chose three different conditions for the clouds at impact, and for each investigated 
adiabatic and radiative flows, giving a total set of six cases. The simplest initial condition involves 
uniform, pressure-bound clouds, immediately adjacent and placed in motion at the start of the 
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simulation in an initially uniform background. Thus, these collisions take place before the clouds 
have formed any structure due to their motions; that is, they are “non-evolved”. Although this 
case is not very realistic, it most closely resembles previous work {e.g., , Stone 1970a, 1970b, 
LMPS, KMW) and, since it produces clean demonstrations of the four stages of CCs, it is very 
useful as a “standard model”. Next, to add some realism and to include nonuniform structures 
in a natural way, we also considered collisions between clouds that had evolved independently 
before impact. Those clouds were, otherwise, identical to the ones used in the “standard”, 
“non-evolved” collisions. We considered two varieties of such “evolved” clouds. First we followed 
collisions between two clouds after each had undergone an identical evolution time (being a 
so-called “crushing time”, defined in equation 2-4), so that the impacting clouds still had a 


mirror symmetry. We alternatively allowed collisions between two clouds of somewhat different 
evolutionary ages. Since during their evolution the clouds become increasingly irregular due to 
KHIs and RTIs, these last collisions involve “non-symmetric” clouds. Thus, we are able in a 
simple, but natural way to begin exploring the role of asymmetry in collisions. Properties of the 
clouds used in each simulation are summarized in Table |^. 

The natural timescale for individual supersonic cloud evolution is the so-called “crushing 
time”, Ter, {e-.g., KMW, Jones et al. 1994) defined as 






Me. 


= 1.3 




pc 
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lOkm s 


-1 


-1 


Myr. 


(2-4) 


During an interval Ter a cloud moving through an external medium develops a bow shock and 
is maximally compressed by an internal shock originating at the front part. An extended, 
low-pressure wake develops behind the cloud. Also on this timescale, KHIs and RTIs start to 
disrupt the cloud as it begins to become decelerated with respect to its background {e.g., Vietri 
et al. 1997). In the absence of magnetic fields, simulations show that clouds disrupt because of 
instabilities on timescales t > Ter- Detailed discussions of the physics of individual, supersonic 
cloud evolution may be found in Doroshkevich & Zel’dovich (1981), Jones et al. (1994), Schiano et 
al. (1995), Jones et al. (1996), Vietri et al. (1997). For the clouds considered here, ~ 5.3 x lO^yr 
in the adiabatic case {Re = 0.4pc) and Ter ^ 2 x 10®y in the radiative case {Re = 1.5pc). These 
are of the same order as the mean time for a cloud to have a collision in the ISM (Stone 1970, 
Spitzer 1978, Hausman 1981). This further supports our feeling that cloud evolution prior to the 
collision must be considered. In our study we begin the individual cloud evolution at t = —Ter 
for “symmetric, evolved” CCs. For “non-symmetric, evolved” CCs one of the clouds begins its 
evolution at t = —l.bTer before the collision event. For these purposes, we define the collision to 
begin {i.e., t = 0) when the bow shocks of the two clouds osculate. 


3. Numerical Setup 






3.1. The Code 


We simulate the CC problem using a 2D Eulerian hydrodynamical code on a Cartesian grid. 
The code we used is based on an explicit “TVD”, conservative finite difference scheme, second 
order accurate both in time and space (Harten 1983, Ryu et al. 1993). Multidimensional flows 
are handled by the Strang-type dimensional splitting (Strang 1968). We accounted for radiative 
cooling in each time step by explicitly correcting the total energy after updating hydrodynamical 
quantities. 

Radiative losses are described generally by the following equation: 

ft = 

where e is the internal energy and C = n^A(n, T) — nE accounts for the net loss to radiation 
against nonadiabatic heating. The cooling term C defines the cooling timescale 

Trad = (3-2) 

When radiative losses are very high, as they can be during the compression phase, the cooling 
timescale, Tradt is comparable to or less than the dynamical time scale which, by the Courant 
condition, ordinarily determines the computational time step. In this case the cooling term is 
labeled “stiff”. There are several ways to handle stiff cooling terms (see, e.g., LeVeque 1997). 
One way is to choose the shorter of Trad and the time step imposed by the Courant condition as 
the computational time step. However, in some situations this choice may lead to uncomfortably 
short computational time steps. On the other hand we could employ Strang’s operator splitting 
approach, where the cooling is computed by multiple steps with its own time step, during one 
computational time step determined by the Courant condition. We chose a third approach in 
which eq. |3-1| is rewritten as 

^ = = (3-3) 

dt 6 '^rad 

This leads to the solution between time steps j and j + 1, 

(3-4) 

For this method the cooling is computed in a single step even though the cooling time scale is 
smaller than the dynamical time scale. The method ensures that the internal energy is always 
positive with the computational time step, At, determined by the Courant condition. If one uses 
the initial value for the radiative correction is, however, only first order accurate in time. 

The radiative cooling function we used, A(n,T), includes free-free emission, recombination lines as 
well as collisional excitation lines with a standard solar metallicity (Z= Zq), whereas the heating 
(r) is provided through ionization and photoionization processes and by cosmic rays (for a full 
description see Ferrara &: Field 1994). We have neglected the effects due to dust grains considered 
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by RFM; this accounts for the slight difference in the ISM phase properties of our model with 
respect to theirs. 

To allow tracking of material that was initially identified with each of the clouds we introduced 
a passive tracer, S, usually referred to as “color” (Xu & Stone, 1995). This quantity (actually, 
one for each cloud) is initially set to unity inside each cloud and zero elsewhere. It represents the 
fraction of material inside each cell that was originally part of one of the clouds. The evolution of 
the color is followed with van Leer’s second-order advection scheme (van Leer 1976). 


3.2. Grid, Boundary Conditions and Tests 

In each simulation the clouds are centered on the x-axis with reflection symmetry assumed 
across this axis. Only the plane y > 0 is included in the computational box. Tests with this code 
show that this more economical grid gives results equivalent to those obtained with a full plane. 
The length scale is adjusted for each case so that Rc = 1.0. In those units, the computational 
domain is x=[0,20] and y=[0,10]. The grid is Cartesian, so that our clouds are actually cylinders. 
The left, top and right boundaries are open. Reflections at these open boundaries are too small 
and too far away from the collision to affect the structure of the flow. Test calculations with 
a computational box twice as large showed no relevant differences from the results we describe 
below. We have explored a range of numerical resolutions, although only the computations at the 
highest resolution are presented here. These involve a 1024 x 512 grid, which provides a resolution 
of 50 zones across the initial cloud radius. We have also compared our results on a uniform 
512 X 256 grid with the Adaptive Mesh Refinement results of KMW. Our simulations are purely 
hydrodynamical, while they added, for convenience a magnetic pressure (but no other magnetic 
effects) to limit compression behind radiative shocks. Still, for the same parameters, we obtain 
results consistent with theirs in terms of main hydrodynamical feature development. 


3.3. Numerical Values, Physical Parameters and Initial Conditions 

The parameters involved in the problem are quite numerous and a comprehensive study 
is beyond our present scope. However, since we are interested in assessing the importance of 
radiative losses, pre-evolution of clouds through the intercloud medium and symmetry of the 
problem, we have decided to restrict as much as possible the parameters’ space by developing an 
accurate model for the multi-phase structure of the ISM, in agreement with observational data. 
We then adopt the most typical values for cloud and WIM physical properties as derived from 
such a model, hoping that they are truly representative of the ISM conditions. In the following we 
give an outline of model assumption. 

As mentioned in §2 we consider clouds that are initially uniform and in pressure equilibrium 
with the inter-cloud medium. We set the initial density contrast y = nc/nj=100, where ric and n* 
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are the cloud and inter-cloud number densities respectively, and ric = 22cm~^ n* = 0 . 22 cm“^); 
the cloud temperature is Tc = lAK (=^ Ti = TdOO-fC). This particular choice is dictated by the 
radiative cooling function adopted in our calculation, by pressure equilibrium assumption, and 
by the density contrast between the two different phases. The equilibrium thermal pressure 
for the ISM turns out to be peq/ks = 1628iir cm~^. Each cloud has an initial Mach number 
M = Vcjcsi = 1.5, where Cgi is the sound speed in the inter-cloud medium. The adiabatic index 
7 = 5/3 (p = [7 — l]e) is assumed throughout the calculations. For the adiabatic cases we set 
Rc = 0.4pc, whereas for the radiative ones, Rc = 1.5pc. With this choice of the parameters we 
have Csi ~ lO/cm and Vc. ~ 15km TcoU ~ 2.6 x IQ^yr for the adiabatic cases. For the 
radiative cases TcoU ~ 9.7 x lO^yr and the radiative cooling time inside the clouds turns out to be 
Trad^ 3-7 X 10^ yr, yielding y ~ 0.38. Finally, the Jeans length associated with the initial clouds 
is Xj ~ 29pc 3> Rc- In the radiative symmetric collisions, the large density increase produced 
during the compression phase, causes a significant reduction of Xj, which becomes comparable to, 
yet still larger than, the vertical size of the clouds. For this reason we have neglected self-gravity 
throughout our calculations (see also KMW). However in a more refined calculation, which 
would take into account other processes like chemical reactions or recombination processes, larger 
compression might be allowed making, as a result, the colliding clouds gravitationally unstable. 

As explained in §2 we have allowed the described adiabatic and radiative clouds to collide 
under three different circumstances. For Cases 1 and 2 in Table |T], uniform clouds are placed on 
the grid in such a way that their initial boundaries are only 2 zones apart at t = 0.0. These are 
the so-called “non-evolved” cases. For Cases 3 and 4 each cloud begins on independent evolution 
at t = —Ter (as appropriately determined by their properties listed in Table H). In these cases 
t = 0.0 is defined as the moment when the bow shocks of the clouds osculate. For Cases 5 and 6 
one of the clouds begins its independent evolution earlier, at t = —I-Btcv Again, the two cloud 
bow shocks come together at t = 0.0. Table lists these details. Animations of each simulation 
have been posted on the World Wide Web site at the University of Minnesota. 

4. Results 

4.1. Collision of non-evolved clouds 

Fig. H shows the four phases defined in §2 for an adiabatic collision (Case 1). At the earliest 
time shown (Fig. la), t = 1.5 TcoU, the collision is near the end of the compression phase. At the 
very beginning of this phase a one-dimensional analysis in the limit of strong shock can still be 
applied to the shocks propagating through the clouds. Theory predicts p ~ Ape and p m)pcvi 
in very good agreement with the numerical values found (Fig. |^. The high pressure in the 
interaction region limits the compression and leads to a fast re-expansion. The compression phase 
lasts longer near the cloud centers, because the cloud column density is maximum along y = 0, 
and because some gas is vertically squirted out from side edges of the interaction region, right 
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after the beginning of the collision (Fig. ||a). The ejected material propagates through the lower 
density WIM, and later on develops features commonly seen in astrophysical jets. In particular 
this slab-jet structure shows a working surface bounded by a shock, a cocoon surrounding the jet 
and apparent KHIs. Nevertheless, since this gas represents a small fraction of the total mass of 
the clouds, it does not affect the development of the collision very much. At the end of this phase 
the gas is highly compressed into a layer much thinner than the initial cloud size (Fig. ||a andH). 

After t = 1.5 Tcozz, as already pointed out in §||, the shocks generated within the clouds 
during the collision enter the WIM and allow the clouds to re-expand. Re-expansion takes place 
supersonically, generating a shock that, with the jet-shock, develops a nearly circularly expanding 
shock-structure on the x-y plane. Inside this structure a reverse shock is generated and the 
re-expanding cloud material begins to accumulate in a dense shell with p r\j 30p, (Fig. I). By t = 
8.2 Tcoii (Fig. Sb) a dense layer is well-formed and, despite its expansion, has become the region 
of highest density (/? ~ 10 — 16pi). It has a nearly circular shape, except for distortions by KHIs 
and RTIs, which eventually form long dense {p ~ 11.5/?z) fingers. The shell expands at pretty high 
velocity {v ~ l.lcsi). After about t =12 TcoH the re-expansion of the shell halts, just before the 
reverse shock passes from the shell into the central, low-density cavity. Then the collapse of the 
shell begins. Fig. ||c shows the collapse phase at t = 37.5 TcoZZ- Eventually, the reverse shock 
reflects off the x-axis and rebounds. But, it is then largely disrupted by refraction in the irregular 
density structure of the collapsing cloud material. Large-scale vortices are generated that through 
KHIs and RTIs hasten the formation of complex filaments evident by t 67.5 Tcoii (Fig. |l|d). 
At the end of the simulation what remains of the two clouds is a very low density central region 
with p mainly between 2pi and 3.5/9* (Fig. |^, surrounded by a complex of filaments, which are 
mixing the original cloud material with the WIM. We deduce that, in Case 1, the likely fate of the 
clouds is disruption and conversion of cloud material into the WIM phase. 

The analogous radiative Case 2 is shown in Fig. During the compression phase (Fig. ^a) 
a fair fraction of the thermal energy is radiated away (Fig. As a consequence, the density 
reaches very high values {p ~ 10^/?*) and re-expansion is much slower than in Case 1. Since the 
re-expansion is so slow now, the reverse shock promptly penetrates all the way back to the impact 
surface and is reflected outwards again. This sequence, much like what happens inside a young 
supernova remnant {e.g., Dohm-Palmer &: Jones 1996), eliminates the central low pressure region. 
Following this, a significant re-expansion along the initial direction of motion gives back to the 
merged cloud material a typical cloud-like aspect ratio. The structure undergoes some vertical 
expansion too, but the most prominent feature in this direction is a thin jet propagating along 
the symmetry plane of the collision. This is shown in Fig. |. At t = 9 T^ou (Fig. |b) , almost 
all the cloud gas is still in a core with high density {p ~ 580/9*), surrounded by lower density 
material with p ranging between 200 and 450/?*, and expanding at a very low speed. A comparison 
of Fig. P? with Fig. |I|b shows clearly that the size of the re-expanded cloud is much larger in 
the adiabatic case than in the radiative one. As a consequence, in the radiative case the slab-jet 
structure is much more distinct. For Case 2 the simulation ends at t = 37.5 TcoZZ- We deduce that 
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in Case 2 the likely fate of elouds is coalescence. The apparently coalesced clouds have evolved 
into an almost circular object of radius ~ 2Rc with densities ranging between about 20pi in the 
inner part and 70pi at the surface. Its total mass is M^ot — 0.84 x 2 x M^., where is the initial 
cloud mass, showing a high efficiency (84%) for the buildup mass mechanism. The edge of the 
newly formed object is sharply bounded, but shows clear signs of KHIs. However, because of the 
high density of its external layer, KHIs will become effective on a timescale much longer than in 
the adiabatic case. The expansion velocity inside the merged cloud at the end is small; namely 
a few xW~^Csi. The net radiative cooling is positive (which means on balance that the gas is 
losing thermal energy) in the outer, denser part of the merged clouds, and negative (which means 
that the gas is being heated up by the background radiation and cosmic rays) in the inner more 
diffuse region, although in both cases the energy gains or losses are not very significant. There 
is a small outward-facing pressure gradient within the cloud concentration, so eventual collapse 
seems likely. Although coalescence seems likely in the near-term, the final fate of the cloud is 
unclear. KMW suggest that the cloud will expand and contract multiple times until pressure and 
thermal equilibrium are reached, developing filamentary structures along the x-axis in the process. 
It might also be possible that the lower pressure inside the cloud induces contraction followed by 
a sufficient increase in the density to turn the cooling function positive. In that case the extra 
pressure due to in-fall could be radiated away and the original cloud density, for which radiative 
equilibrium holds, might be approached. In addition, because of the highly dense external layer 
around the newly formed structure, we do not expect in the radiative case that RTIs will be as 
disruptive during the collapse phase as for the adiabatic case. 


4.2. Collision of evolved clouds: symmetric cases 

Fig. ^ shows the initial conditions for the adiabatic collision of two evolved clouds (Case 3). 
At t = 0, when the cloud bow shocks just touch and each cloud has evolved through one crushing 
time, the density changes smoothly through the clouds, ranging from 10/3j, at the back, to 150 pt 
at the front part. The x-component of the velocity follows the same pattern, being higher at the 
front of the cloud than at the rear, although the range of this variable is much smaller {e.g., Jones 
et al. 1994). On average the clouds have an individual speed corresponding to a Mach number, 
M = 1.35, relative to the WIM. There are a number of differences introduced into the interaction 
between the clouds by allowing for prior evolution. The most important are the presence of 
bow shocks and incoming gas motions within each cloud wake once the clouds collide. After the 
approaching clouds encounter each other’s bow shock, reverse shocks (which act as secondary bow 
shocks) are generated. In the adiabatic interaction (Case 3) these shocks substantially affect the 
clouds. In fact at t = 2.2 TcoIU right before the cloud bodies encounter the bow shocks, their 
x-width I ~ 1.5i?c and pmax ~ 150, but at t = 3 TcoU, after the bow shock-cloud collision begins, 
p ~ 250pi ( Pmax^ 260pj) and I ~ 1.2i?c- There is further compression so that at t ~ 3.8 TcoIU 
right before the cloud bodies impact each other, p ~ 380pi in the compressed front layer, and 
£ ^ Rc (Fig. ^). So, the maximum compression reached during the collision ( t ~ 4.5 Tcou) 
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( Pmax^ 1200Pi) is much higher in Case 3 than in Case 1, although the pressure enhancement 
is about the same, and both collisions are adiabatic. If the clouds were self-gravitating, such 
differences might become important to the possibility for triggering star formation out of such 
collisions. On the other hand, the jet-like and the thin shell structures develop in roughly the 
same pattern as for Case 1. By comparison to the non-evolved collision, however, more cloud 
material remains in a core structure, almost to the re-expansion phase. This is evident in Fig. lb, 
at t = 11.2 Tcoii, when contrasted with Fig. |^. The same comparison also shows that the shell 
structure in Case 3 is more irregular than in Case 1. This is in fact due to the lower density of 
the shell of Case 3, which allows a quicker development of RTIs and KHIs. By the time shown in 
Fig. Ic expansion along the x-axis has been reduced compared to the y-direction. That results 
from the interaction between the expanding cloud material and inflowing WIM within the wakes 
of the two clouds. The ram pressure of the wake flows is strong enough to affect significantly 
the expansion along this path. Indeed at t ~15 TcoU a standing, reverse shock is well formed 
and deflects the expanding cloud material away from the x-axis. This in turn becomes strongly 
sheared and filamentary due to KHIs. A substantial fraction of the cloud material is ejected 
so that it cannot join the collapse. The shock structures in this case are very complex, since 
they involve interactions with the pre-existing bow shocks and tail shocks and a generally more 
complex density and pressure structure as the collision begins. Nevertheless, as our simulation of 
Case 3 ends {t = 22.5 Tcou), the dominant, irregular shell structure resembles qualitatively that in 
Case 1. Subsequent evolution in Case 3 should follow a pattern similar to that in the analogous 
non-evolved Case 1. In particular we expect disruption of the clouds in both cases. 

In the radiative evolved case (Case 4) the individual clouds show a similar qualitative 
structure as in Case 3, but now their compression due to initial motion through the inter-cloud 
medium is much higher ( Pmax^ 14:00pi). The additional compression results from enhanced 
radiative cooling induced by the shock compression within the cloud. We note that the bow-shock 
compression phase that was important to Case 3 turns out not to produce a very significant effect 
for Case 4. That is because the compression brought on by radiative losses prior to the encounter 
is already very large. In addition, since the speed at which the bow shock penetrates each cloud 
scales inversely with the square root of density, in this case the bow shocks barely penetrate into 
the clouds bodies before they collide, producing only little pre-compression with respect to the 
adiabatic case. As a result, the compression reached during the collision is only slightly higher ( 
Pmax^ 4 X 10^Pi) than in the non-evolved case (Case 2). The re-expansion phase in Case 4 follows 
pretty much the same pattern as in Case 2, except that now it is substantially slowed down in 
some directions by the action of the wakes behind each cloud, as we noted also for Case 3. In fact, 
as shown in Fig. by the end of this simulation ( t = 22.5 Tcou) the clouds have merged into a 
dense structure of size ~ 2Rc, which is 2/3 of the size of the merged core in Case 2. Expansion in 
the y-direction is relatively free and leads to a KH unstable jet as in the other symmetric collisions. 
However, re-expansion along the x-axis is strongly inhibited by inflowing wake material, so that 
the velocity of expansion along the x-axis is only about 1/3 as large, compared to Case 1. As a 
result the expanding gas collects in a high density cloud “rim” (p ~ 400pj, Pmax^ 440pj). The 
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wake material forms a standing, outward facing, “accretion shock” outside the cloud structure. 
Examination of the color variable shows a very good correspondence between the high density 
material visible in Fig. and cloud material. Thus, the figure shows that wake material, once it 
impinges on the cloud, joins the outflow of cloud material in the y-direction. In fact, it appears 
that the inflowing wake material and associated shocks are responsible for driving the vertical 
outflow of cloud material along the jet and for producing the KHIs that have generated the large 
eddies evident in Fig. |^. In the inner cloud the density ranges between 100 and 200 pi, which is 
higher by almost a factor 2 with respect to the non-evolved case. Even though there is a weak 
pressure gradient pushing vertically in the central condensation, it seems very likely that in an 
extension of this simulation the main core of the merged clouds would remain intact. So, we judge 
these clouds to he coalesced. 


4.3. Collision of evolved clouds: asymmetric cases 

As mentioned earlier we chose a simple, but natural way to explore symmetry breaking in 
the collisions just discussed. That is, we collided clouds that were identical when set into motion, 
but which were differently evolved at impact. On a timescale ~ Ter the compression substantially 
deforms the clouds, while KHIs and RTIs will produce irregular cloud boundaries. Since those 
features are highly time dependent, two clouds of even slightly different dynamical ages will lack 
mirror symmetry. In our asymmetric simulations one cloud (Cl hereafter, and on the left) was 
aged by 1 Ter and the other (C2 hereafter, and on the right) by 1.5 as their bow shocks came 
in contact at t = 0.0. The older cloud had in general a denser front part and a smaller velocity at 
impact (Jones et al. 1994). However the aspect ratio (length to height ratio) developed by the 
cloud during its motion through the WIM is strongly related to the adiabaticity of the gas. Indeed 
it increases in the adiabatic case, but decreases in the radiative one, so that the radiative cloud 
grows denser and more compact as it evolves (Vietri et al. 1997). This turns out to have a major 
impact on the survival of clouds in asymmetric CCs. As in the previous evolved cases, the two 
clouds undergo bow shock-compression before colliding bodily. In the asymmetric interactions 
the clouds have different speeds and shapes, and are located at different distances from their bow 
shocks. As a result they no longer experience compression simultaneously. In addition, since C2’s 
bow shock is weaker than Cl’s, bow-shock compression for C2 turns out to be stronger than for 
Cl. In the adiabatic Case 5, right before the direct collision. Cl has p ~ SOOpi and £ ~ i?c and C2 
has p ~ 370/9j and ^ ~ 0.6i2c- The compression phase is shorter than in Cases 1-4, and not all 
of the kinetic energy is converted into thermal energy (Fig. ^). After the collision, the younger, 
more compact cloud Cl, maintains its identity longer than C2. As shown in Fig. ^ at t = 8.2 
TeoU Cl still has a dense core with p ~ 70/?j surrounded by a layer with p ~ 40/?i. By contrast, C2 
is being stretched and torn apart and, soon after, is mostly converted into the WIM. Nevertheless 
Cl is undergoing rapid re-expansion too and, at t = 13.5 Teoiu the density is lower than 15 pi 
everywhere. Although it is much more irregular than Cases 1 or 3, in Case 5 we can still recognize 
a clear pattern in the evolution of the re-expansion phase, with the formation of an expanding 
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shock wave along with a reverse shock and a low density shell, heavily affected by KHIs and 
RTIs. So, at the end of the second, re-expansion phase, part of the cloud gas has already been 
converted into WIM, whereas the remaining part is forming an irregular hlamentary structure with 
p ranging between 5 and lO/Oj. As in the symmetric-evolved collisions, the interaction between 
re-expanding cloud material and inflowing wake gas significantly influences the expansion. As for 
Case 3 the wake confines expansion along the x-axis through a standing, inward-facing shock. 
Also, as in Case 3, the expanding cloud material is deflected around the wake in a strongly sheared 
manner. The broken symmetry in Case 5 allows C2 material to expand much more easily, however, 
because its greater y-extent at impact effectively “launches” it over the wake of Cl. In addition, 
since C2 has a lower column density just above the x-axis, it recoils in response to the collision. 
Thus, a substantial fraction of C2 passes over its own wake after the collision. Those effects add 
considerably to the disruption of cloud material into hlaments and probably hastens its dispersal 
into the WIM. The resulting flow is very complex and highly vortical. So, we expect the collapse 
phase to be very ineffective at collecting together material into new clouds. Thus, this case seems 
clearly much more disruptive that any of the previous ones. 

On the other hand, the radiative Case 6 is characterized by completely new features. At t 
= 2.2 Tcou the region between the two shocks formed by the collision, appears to be strongly 
distorted, resembling the structure that develops during the “bending mode instability” (KMW, 
Vishniac 1994, Hunter et al. 1986). As shown in Fig. |^a, at t = 3 TcoU the layer between the 
two shocks has grown more corrugated and the front parts of the clouds are following the same 
pattern. Although the bending mode is unstable, the disruption of the impact surface is primarily 
due to the gross irregularities of clouds’ shape, which dominate the determination of the evolution 
of the structure of the flow. Fig. ^a shows also that the upper part of the left (Cl) cloud is 
about to break off and pass over the other cloud (C2), carrying some of C2’s gas with it. As 
already pointed out, unlike the adiabatic case, in the radiative case the older cloud C2 has a more 
compact structure, which makes it more solid and more resistant to the collision than Cl. At t 
= 9 Tcoii (Fig. 0b), the material from Cl passing over the top of C2 has expanded again into a 
distinct cloudlet, with density ranging from p ~ lOOpi at the front, to p ~ 5pi at the tail, and 
with still slightly supersonic velocity {vx > Csi). Subsequently, that cloudlet becomes strongly 
decelerated and suffers KHI- and RTF induced destruction, as seen previously for individual 
clouds {e.g., Jones et al. 1994). Although the cloudlet leaves our grid before it is destroyed, it 
seems fairly clear that it will dissolve into the WIM. The wake of the cloudlet is still visible on 
the far right side in Fig. |^. On the other hand the core of C2, which by t = 3 TcoU has a very 
long tail and a dense, but distorted front, passes through the remaining part of Cl, emerging after 
t Ri 7.5 Tcoii with the main body of C2 accreted. This outcome resembles those one would expect 
in collisions of two clouds of quite different sizes and densities. It shows that the smaller and 
denser cloud is able to pass through the larger and more diffuse one, sweeping its gas and finally 
breaking it up into two major pieces (Gilden 1984, Kimura &: Tosa 1996). In these calculations we 
show that the same fate can actually occur when two clouds with approximately the same initial 
characteristics but with slightly different morphologies collide. 
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As in the other collisions between clouds that are followed by wakes, an “accretion” shock 
forms to the right of C2. However, it is driven off the grid to the right before the simulation ends. 
The interaction is dominated in this case by the high concentration of material in Cl. 

The merged remnant of Cl and C2 that was formed by t =7.5 TcoU has by the end of 
the simulation, t = 30 TcoU, evolved into non uniform filamentary structures, characterized by 
irregular motion (Fig. |^c). Very little mixing between the two original clouds has actually taken 
place; rather one has passed through the other. On the other hand, considerable entrainment of 
WIM gas has taken place through the action of eddies generated during the collision. The higher 
density features apparent on the left and top perimeter of the main cloud visible in Fig. ^ are, in 
fact, the remnant merged cores of the original clouds, whereas the rest of the main cloud at this 
late time contains a strong mix of entrained material, or has been ejected from the grid. At this 
time the main cloud is being slowly stretched along the x-axis (vx ~ O.lc^i at the edges), while 
large eddies on the top and downward pointing pressure forces are being effective at reducing 
its height. It seems likely that the single dense region visible in Fig. will be bisected into 
two before too long. The final outcome may be two distinct clouds formed largely from material 
originally in C2, which was the more compact of the original pair. Thus the outcome is completely 
different from that of Cases 2 and 4 and we conclude that the two clouds are destroyed by the 
collision and converted into several filamentary structures. 


5. Discussion 

We have investigated for the six cases summarized in Table |I] the collision of diffuse HI clouds 
in a multi-phase medium. Our objective is to understand such issues as the likely fate of clouds 
after collisions, including the conditions for coalescence and the fraction of the initial kinetic 
energy radiated away. In the previous section we outlined the basic dynamical evolution of each 
collision and the ultimate fate of the clouds. It was clear from those examples that the fate of 
colliding clouds depends strongly on the symmetry of the interaction and also on the degree to 
which the initial kinetic energy is radiated away before the clouds begin their re-expansion. To 
clarify and expand on those issues we now review the main points, separately for the adiabatic 
and radiative cases. 


5.1. Adiabatic Cases 

Adiabatic collisions generally appear to result in cloud disruption, with most of the gas 
converted into the WIM phase. This point is made clearer in Fig. which shows a plot of various 
properties characterizing individual clouds involved in adiabatic collisions. In particular, using the 
color variables, we can follow the kinetic and thermal energies of each cloud, normalized to the 
initial cloud total energy, as well as the center-of-mass coordinates Xcm and ycm for each cloud 
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relative to the point of first contact (see Jones et al. 1996 for a mathematical definition of Xcm and 
ycm)- From a comparison of panels (a) and (b) in Fig. we see that in all adiabatic Cases 1, 3 and 
5 the kinetic energy of the clouds Flfcm = is initially converted mostly into thermal energy 

Ether = pyHi — 1) during the compression phase, through the action of the main, outward moving 
shocks as they pass through the cloud bodies. That stage is immediate for the “nonevolved” 
collision, but is delayed, of course, for the “evolved” cases, since t = 0 corresponds to the moment 
when the bow shocks touch, rather than when the clouds first touch. During the re-expansion 
phase, some of this thermal energy is converted back into kinetic form, since the rapid expansion 
of the “blast wave” into the WIM, reduces the pressure around the merged clouds. But, because 
the cloud shocks have generated entropy, the cloud gas has substantially more thermal energy at 
the end of each simulation than at the start, despite considerable volume expansion. This effect 
enhances the tendency for the cloud material to be converted into the WIM phase following the 
collision. Also, these plots show that the normalized kinetic energy decreases gradually with time; 
this is due to the irreversible work done by the expanding gas on the surrounding background. A 
similar trend was observed in the 3D simulations of two colliding gas streams by Lee et al. (1996). 
It is important to notice that at late times {t >8 Tcou) some energy decrease is also due to escape 
of matter from the computational domain, particularly the top. The center-of-mass positions 
provide a quantitative measure of this expansion as shown in panels (c) and (d). In particular, 
note for Case 1 that the maximum in 7/cm around t >8 TcoU corresponds to the time when the 
dense shell reaches the top of the grid. After that time 7/cm represents only mass the remaining 
within the grid. Panel (c) also illustrates quantitatively how cloud material in the symmetric 
collisions is more efficiently removed from the central interaction region than in the asymmetric 
case. Notice for Case 5 that some cloud material remained in the interaction region interior. In 
panel (d) we also see the clear difference in y-expansion of the individual clouds in the asymmetric 
Case 5. The more evolved cloud, C2, which begins the simulation having a greater height and 
lower column density along central impact axis, is obviously disrupted by this measure. The more 
compact cloud. Cl, remains compact in this dimension, on the other hand. 


5.2. Radiative cases 

In the radiative cases, as long as the initial geometry is symmetric, the two colliding clouds 
merge and it seems likely that the conditions for a new massive cloud to form do exist. As shown 
in Fig. 1^, the kinetic energy of the two clouds is converted into thermal energy and soon radiated 
away. In deed, unlike the adiabatic cases, after the compression phase there is no re-enhancement 
of the kinetic energy, as is shown by a comparison of panels (c) and (d) in Fig. ^ and p. As a 
result the re-expansion of the gas is much reduced with respect to the adiabatic case. At the end 
of the evolution the thermal energy of the clouds has increased, due to heating processes occurring 
during the re-expansion phase (Fig. P). Also, in the evolved (symmetric) cases, the density at the 
compressed layer formed during the collision is higher than in the non-evolved case and this could 
be important for triggering star formation processes. Furthermore, the wakes developed during 
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the cloud evolution limit the re-expansion, so that after the collision in Case 2 the merged clouds 
end up forming a smaller and denser region. In cases where self-gravity is dynamically important, 
the wake effect could play a fundamental role, reducing the expansion velocity of the material 
below the escape threshold, thus making it possible for the two clouds to build up a new larger, 
gravitationally bound structure (see, also, Dinge (1997) for similar behavior in a single, moving, 
self-gravitating cloud). 

The outcome of CCs is very different for the asymmetric, i.e. differently evolved clouds, even 
though our collisions all involve clouds of equal mass and head-on collisions. In this case we obtain 
fragmentation of the “younger” cloud and the formation of a small expanding cloudlet eventually 
dissolving into the WIM. As shown in Fig. |^, in this case only part of the initial kinetic energy 
is radiated away during the collision. As a result, the remainder of the merged cloud material 
re-expands and, through interaction with the wake flow, becomes concentrated into dense clumps. 
At the end most of the cloud material has been converted into WIM, its directed kinetic energy 
partly being lost to radiation, partly to turbulent gas motion. KMW have shown that by slightly 
perturbing the surface of one of the two colliding clouds, the bending mode instability causes 
fragmentation rather than coalescence. In our study the asymmetry in the problem is due to a 
slight difference in the evolutionary ages of the clouds (0.5 Ter)- As already pointed out, its effect is 
comparable to that of an off-center collision. That is, the mass and momentum of the two clouds at 
impact have significantly different y distributions. Those structural differences of the clouds make 
the outcome of the collision very different from the symmetric cases. From this fact we conclude 
that supersonic radiative collisions of clouds with the same gross characteristics (mass, speed and 
structure), but without high symmetry, are disruptive and generate irregular filamentary clumps. 
On the other hand previous results show that collisions of clouds with major structure differences 
(density, size), are likely to be disruptive as well (Hausman 1981, Gilden 1984, Kimura &: Tosa 
1996, RFM). Our results support this. In addition, our results show that the adiabatic case, which 
applies to small clouds, will be highly disruptive, even if the clouds are identical and the collision 
head-on. Also, as pointed out below, off-center collisions, which are the most common, certainly 
are very unlikely to produce coalescence of clouds, at least for non-self-gravitating objects. 

These results have an important impact on the ISM of galaxies. In general, supersonic, 
gasdynamical CCs tend to be disruptive, so that any model to explain cloud evolution and mass 
spectrum has to take these findings into account. If colliding clouds are sufficiently electrically 
conductive, then the presence of a large-scale magnetic field may play an important dynamical 
role in the interaction. It is not obvious, however, that our main conclusion will be altered. We 
have fully MHD simulations underway that will address that point in a separate report. 

It is worthwhile to mention that our asymmetric calculations can provide some insights for 
off-center collisions as well, at least for those with small impact parameter (6 <C Rc)- In fact, 
with regard to survival of the clouds, the most important implication for small impact parameter 
off-center collisions is probably asymmetry. For these cases the results of our asymmetric 
calculations may apply as a guide, at least for clouds of comparable mass. However, when the 
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impact parameter b ^ Rc, only part of the cloud is involved in the collision and our calculations 
are not appropriate anymore. LMPS have investigated off-center, isothermal collisions for very 
massive clouds, including gravity in their calculations. They show that, for high relative velocity, 
the colliding parts of the clouds soon re-expand and disperse after the compression phase, whereas 
the outer parts (which do not get involved in the collision) proceed unhindered and form two new 
small clouds. On the other hand, for low relative velocity, the colliding clouds coalesce, whereas 
the outer gas motion is deflected into a circular pattern. As a result rotating bound systems 
and bars form. For smaller non-self-gravitating clouds, we think that, in addition to the relative 
velocity, the adiabaticity is a key parameter. If b < Rc the collision can probably be classified in 
similar terms to a highly asymmetric one. On the other hand, if 6 > Rc, the LMPS calculations 
suggested that the outer part of the clouds are torn apart during the collision. For an adiabatic 
collision, pressure waves in the remaining clouds generated during the encounter might be able to 
make the clouds expand and disperse into the WIM. However, in a strongly radiative encounter 
these waves could be damped away by radiation and the cloud cores might survive, although the 
clouds themselves would turn out quite distorted Finally, if 6 ~ 2Rc, the collision will produce 
only minor perturbations on the cloud structures. 

Real clouds are, of course, three dimensional. CC three-dimensional calculations have been 
performed by Hausman (1981), LMPS and Lattanzio & Henriksen (1988, LH hereinafter); these 
authors, using a smoothed particle hydrodynamics code, investigate on the effects of several 
parameters (relative velocity, cloud mass ratio, impact parameter and so on) on coalescence. 
Hausman’s calculations are strongly affected by the limitations of his computational means. 

In particular he used a resolution as low as 100 particles per cloud, and his calculations show 
unphysical particle interpenetration. For this reason, for example, in his run 1, Hausman finds 
a ratio Pmax/Pext of only 3.55, much less than 50 as found by LMPS for the same case. As a 
result, little conversion of kinetic energy into thermal energy takes place; it is probably for this 
reason that, in most of his runs, Hausman finds that the cooling is efficient enough to ensure 
isothermality (Hausman 1981). Since the set of cases studied by Hausman is very similar to that 
presented by LMPS, we will neglect to go in more details about his results. LMPS employed 
a better version of smoothed particle hydrodynamics code (Monaghan & Lattanzio 1985), and 
higher resolution (~ 2000 particles per cloud). We have already summarized their results regarding 
off-center collisions. They also found that symmetric head-on collisions generate a single rapidly 
re-expanding or collapsing cloud depending on whether the initial clouds are gravitationally stable 
or unstable respectively. In asymmetric (but still head-on) cases however, they concluded that 
even when a cloud is marginally stable to gravitational collapse, the collision with a smaller 
cloud (ratio of masses larger than 2.5; their clouds had the same density), does not initiate the 
instability. Finally LH further investigated this problem showing how spin and orbital angular 
momentum are important in determining the outcome of CCs. The accuracy of these results, 
although useful and to some extent in agreement with previously cited works and our own, is 
still limited by low resolution (~ 2000, 3000 particles) and some assumptions which are not 
always appropriate. Some of these were mentioned by the authors themselves. In addition, based 
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on Hausman (1981) results, both LMPS abd LH assumed isothermal clouds, using 7 = 1; but 
as already pointed out, this approach is not correct because of the fundamental dynamical role 
played in the re-expansion phase by thermal energy stored during the compression phase. That 
aspect is apparent in our results. In our radiative head-on symmetric simulations, which allow for 
a release of energy through radiative processes, the collisions produce a new merged stable cloud 
instead of a rapidly re-expanding cloud as found by LMPS. Also LMPS and LH do not include 
in their calculations an intercloud medium. However, as we have noted, the interaction of cloud 
material with the intercloud gas, particularly with shocks and cloud wakes, significantly affects 
the evolution of the collision. 

Useful insights about how an additional degree of freedom in this problem might modify 2D 
results may be provided by 3D studies of some related problems. For example, Xu & Stone (1995) 
examined the evolution in 3D of a gas cloud overrun by a plane shock. They found behaviors 
qualitatively consistent with 2D simulations of that problem {e.g., Bedogni & Woodward 1990; 
Jones &: Kang 1993; KMW). In many respects a shocked cloud is similar to a supersonic cloud in 
its evolution (Jones et al. 1994). Lee et al. (1996) have carried out a 3D study of two colliding 
gas streams. That is similar in some ways to collisions between clouds. They found, as we do 
for adiabatic collisions, that the bulk of the kinetic energy is converted into thermal pressure, 
and that this causes the colliding material to expand as it drives a shock into the ambient 
medium. Generally, the extra dimension will allow more complex motions to develop, and some 
considerations will depend quantitatively on the third dimension, but we see no evidence in the 
existing literature that the general conclusions of the present work will be invalidated when it is 
included. 


6. Conclusions & Summary 

To summarize, in this paper we have found the following results for the gasdynamical 
collisions between two mildly supersonic interstellar clouds, M > 1.5, of comparable mass: 

1. Supersonic CCs are most often disruptive. In particular adiabatic collisions, which involve 
small clouds {Rc < OApc, for the standard WIM parameters assumed in our calculations), 
turn out always to be disruptive. 

2. For completely symmetric collisions, strong radiative energy losses can, however, lead to 
coalescence of the two clouds. In fact emission of radiation reduces the thermal energy 
stored during the compression phase, preventing a vigorous pressure driven re-expansion. 

3. Asymmetry in the clouds at impact greatly enhances the tendency for clouds to be disrupted 
during the interaction, even when radiative cooling is strong. This is true even for a very 
modest asymmetry. In the adiabatic, asymmetric case the clouds are almost immediately 
dispersed in the WIM. In the radiative case new filamentary structures are produced out of 
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the initial cloud material. 

4. Future numerical work should not neglect the importance of allowing the clouds to develop 
a self-consistent structure, especially bow shocks and wakes, since these features strongly 
influence the interaction and add important hydrodynamical features. In particular, bow 
shock interaction of the colliding clouds produces higher compression, particularly in the 
adiabatic case. On the other hand the wakes behind the clouds reduce the re-expansion 
along the x-axis, increasing the probability for coalescence. 
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Table 1. Summary of 2D-HD Cloud Collisions Simulations 


Case 


Rc 

'^coll — Rc / Vc 

'^CI - ‘^Rc -sj^l^c 


Clouds ages® 

End Time® 

1 

adiabatic 

0.4 pc 

0 

X 

5.3x10® 22 

3 

0 Tcr~ 0 Ter 

67.5 Tcoll 

2 

0.38 

1.5 pc 

9.7x10^2/ 

2.Ox 10® 22 

3 

0 Tcr~ 0 Ter 

37.5 Tcoll 

3 

adiabatic 

0.4 pc 

2.6x10"^ 2/ 

5.3x10® 22 

2.7 

1 Ter~ 1 Ter 

22.5 Tcoll 

4 

0.38 

1.5 pc 

9.7x10"^ 2 / 

2.Ox 10® 22 

2.8 

1 Ter~ 1 Ter 

22.5 Tcoll 

5 

adiabatic 

0.4 pc 

2.6x10"^ 2 / 

5.3x10® 22 

2.4 

1 Tcr~ 1-5 Ter 

22.5 Tcoll 

6 

0.38 

1.5 pc 

9.7x10"^ 2 / 

2.Ox 10® 22 

2.7 

1 Tcr~ 1-5 Ter 

30.0 Tcoll 

^All models have used 7 

= 5/3, X = 

= pc! pi = 100, equilibrium pressure Peg /fcs = 

1628 K cm ®. Also, at equilibrium. 

we have Ti 

= 7400K and 

m = 0.22 

cm~^ for the inter-cloud medium and Tc = 

747^ and ric = 22 cm~^ inside the 


clouds. Furthermore, all the computations were carried out on a rectangular grid with size Nx = 2Ny = 1024 
corresponding to a resolution of 50 zones through the cloud radius. Left, top and right boundaries are free whereas 
reflection is assumed at the bottom. 

r/ = Tradj Tcoll 

‘^This is the relative Mach number and is referred to the inter-cloud sound speed, Cai ~ lOfcm 

‘^The left and right hand columns refer to the left and right hand clouds. Cl and C2, respectively. 

®The end time is expressed in terms of collision time TcoU, and represents the total time from the beginning of the 
collision. 
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FIGURE CAPTIONS 


Fig. 1.— Inverted grayscale images of tanh(log(p)) for Case 1; vector display of the velocity field 
superimposed. Panel (a) (top) shows the compression phase at t =1.5 TcoU with outflow at the side 
of the contact discontinuity. Panel (b) (bottom) shows the re-expansion phase at 8.2 TcoIU with the 
formation of a dense shell-like structure. Panels (c) and (d) correspond respectively to the collapse 
phase at t =37.5 TcoU and to the dispersal phase at t =67.5 TcoU- The dramatic development of 
KHIs and RTIs is evident. 

Fig. 2.— Cut through the grid along the x-axis (y=0.3Rc) for Case 1. The top and bottom panel 
show log-plot of density and pressure respectively. Solid lines refer to t =1.5 TcoIU dot lines to t 
=8.2 TcoU, short-dash lines to t =37.5 TcoU and long-dash lines to t =67.5 TcoU- 

Fig. 3.— Same grayscale as in Fig. [I], but for Case 2. (a) (top) The compression phase at t = 
1-5 TcoU, and the ejection of material along the collision plane of the crushed clouds, (b) (center) 
The re-expansion phase at t = 3.7 TcoU', the (slab) jet-like structure is well-formed and is showing 
features common to astrophysical jets, (c) (bottom) The structure at t = 37.5 TcoU, when the two 
clouds have merged into a structure with a dense rim. By this time KHI structures are evident 
along the jet and on the perimeter of the merged cloud. 

Fig. 4.— Same grayscale as in Fig. ||, but now for Case 3. (a) Conditions when the bow shocks 
generated by the clouds osculate, just prior to impact, (b) The two clouds a little before the actual 
cloud-body collision. The clouds have been compressed by the bow shocks, (c) The re-expansion 
phase at t = 11.2 TcoU- The two clouds are still distinguishable; also a low density layer has formed 
around them and is undergoing strong ablation by KHIs. 

Fig. 5.— Same grayscale as in Fig. ||, but now for Case 4, with velocity vectors superimposed. 
This image shows the hnal stage of the re-expansion at t = 22.5 Tcou- The two clouds have merged 
into a new structure with density decreasing from the bottom to the top. Delicate density features 
have been produced by KHIs, enhanced by in flowing wake material deflected along the sides of the 
clouds. 

Fig. 6.— Same grayscale as in Fig. ||, but now for Case 5, with velocity vectors superimposed. 
This image shows the complicated re-expansion phase at t = 8.2 TcoU- The cloud originating on 
the left (Cl) is still recognizable, whereas C2, coming from the right, has been strongly distorted 
from its original form. 

Fig. 7.— Same grayscale as in Fig. ||, but now for Case 6. (a) The compression phase at t = 3.7 
TcoU- The contact surface is quite irregular and asymmetric, foretelling of the disruption to follow, 
(b) The structure at t = 9 TcoU showing the formation of a cloudlet breaking off from Cl. (c) 
Density with superimposed velocity vectors, representing the re-expansion phase at t = 30 TcoU- 
A large irregular structure has formed and eddies are shredding the top part. 
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Fig. 8.— The four panels show plots of the kinetic energy (a - top left), thermal energy (b - bottom 
left), center of mass X-cm (c - top right) and Y-cm coordinates (d - bottom right), as a function 
of time, for each cloud involved in the collision. For Cases 1 and 3 the two colliding clouds are 
identical to each other, so only one is displayed in each case. Solid lines correspond to Case 1, 
dotted lines to Case 3, short-dash lines and dot-long dash respectively to cloud Cl (left) and C2 
(right) in Case 5. The delay of the features of Cases 3 and 5 with respect to the non-evolved Case 
1, is due to the delay in the former cases between impact of bow shocks and cloud bodies. 

Fig. 9.— Same as in Fig. but now solid lines refer to Case 2, dot lines to Case 4, short-dash lines 
and dot-long dash respectively to cloud Cl (left) and C2 (right) in Case 6. Taking care to note 
the significant difference in the Y-axis scale, A comparison with panels (a) and (b) of Fig. shows 
that in the radiative cases much less kinetic energy is converted into thermal form. In addition, as 
shown by the and Ycm panels ((c) and (d) respectively) and the significant difference in the 
y-axis scales of the two figures, re-expansion is much reduced. 
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